Minimal Agent Based Model for Financial Markets 
II: Statistical Properties of the Linear and 
Multiplicative Dynamics 

V. M. Cristelli\ L.Pietronero^'^ A. Zaccaria^ 

^ Universita "La Sapienza", P.le A. Moro 2, 00185, Roma, Italy 
Centro "E. Fermi", Compendio Viminale, 00184, Roma, Italy 
^ ISC-CNR, V. dei Tam-ini 19, 00185, Roma, Italy 

December 5, 2008 



Abstract 

We present a detailed study of the statistical properties of the Agent 
Based Model introduced in paper I [T] and of its generalization to 
the multiplicative dynamics. The aim of the model is to consider the 
minimal elements for the understanding of the origin of the Stylised 
Facts and their Self-Organization. The key elements are fundamental- 
ist agents, chartist agents, herding dynamics and price behavior. The 
first two elements correspond to the competition between stability and 
instability tendencies in the market. The herding behavior governs the 
possibility of the agents to change strategy and it is a crucial element 
of this class of models. The linear approximation permits a simple 
intepretation of the model dynamics and, for many properties, it is 
possible to derive analitical results. The generalized non linear dy- 
namics results to be extremely more sensible to the parameter space 
and much more difficult to analyse and control. The main results for 
the nature and Self-Organization of the Stylised Facts are, however, 
very similar in the two cases. The main peculiarity of the non linear 
dynamics is an enhancement of the fluctuations and a more marked 
evidence of the Stylised Facts. We will also discuss some modifications 
of the model to introduce more realistic elements with respect to the 
real markets. 



1 Introduction 

In the preceding paper (paper I [T]) we have introduced a minimal agent 
based model to discuss the origin and self-organization of the Stylised 
Facts (SF), which are a common characteristics of all price time series. 
The model is based on four essential elements (first introduced by Lux 
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and Marches! (LM) [21 E]): fondamcntalists (F), chartists (C), herd- 
ing, price behavior. Our aim, in the construction of the agent based 
model, has been to look for the maximun simplification and reduction 
of free parameters in order to achieve a detailed understanding of the 
dynamics of the system. This is achieved by a series of simplifications, 
discussed in paper I [T] which, however, maintain the essential elements 
of the model. 

The construction of an agent-based model represents necessarily a com- 
promise between simplicity and realism, many of the models which have 
been proposed, are able indeed to reproduce some of the SF. However, 
often the origin of these properties remains hidden in the many pa- 
rameters of the models. In this perspective the addition of too many 
elements ,which make the model more realistic, may in the end also 
represent an element of confusion for a detailed understanding of the 
dynamics of the model. This point of view has been argued very clearly 
in three recent reviews [H El E] ■ These authors in fact point out that 
it would be highly desirable to introduce a model which is as simple as 
possible but it is still to able to repoduce the SF. This is precisely the 
point of view we have adopted in this and in previous paper where we 
describe in detail this minimal model. Once the essential elements are 
clarified it is certainly possible to generalise the model by introducing 
more realistic elements, and this will be the object of our future work. 
In paper I [1] we have introduced and discussed the essential features 
of the model. The SF are shown to correspond to finite size effect, 
with important conceptual and pratical implications. We propose that 
the self-organised state is linked to a threshold in the agent activities 
which strongly modulates the number of active agents. This leads to a 
feedback mechanism which triggers spontaneously the system towards 
the self-organized intermittency state. This state is not really critical, 
in the sense of statistical physics, because it is related to finite size 
effects which can, however , occur at different scales. 
In the present paper we are going to examine in details the statistical 
properties of our model in various directions. We will try, as much as 
possible, to derive analitical results. This possibility is much easier if 
one considers a linear dynamics for the price instead of the more re- 
alistic multiplicative dynamics. The two cases essentially coincide for 
small price variations and we will consider in details their relation and 
possible discrepancies. 
The paper is organized as it follows. 

Section 2 provides some analytical results in two important regimes of 
the model and in section 3 we investigate the non trivial diffusional 
properties. In section 4 we focus on the behavior of the tails of the 
return probability density function (pdf) when the time lag at which 
the calculation of returns is performed becomes larger and larger. The 
results of Section 4 suggest a simple approximation that explains how 
the tails of return pdf are generated by the model. The approximation 
also allows us to test for which SF property the population dynamics 
plays a crucial role and this will be discussed in Section 5. In section 
6 we introduce an additional RW performed by the fundamental price 
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and this will lead to interesting crossover phenomena. In Section 7 
we consider the multiplicative version of this minimal ABM, and we 
discuss some non trivial aspects with respect to the linear dynamics. 
We conlude in section 8 with a short summary of the results of our 
analysis also with respect to the empirical results, with a discussion 
of future developments of the model. In appendix A we discuss the 
implication of the use of diverse definition of the returns. In appendix 
B we discuss the statistical significance of dimension of a dataset with 
respect to various properties. Finally in appendix C we discuss data 
collapse for the return pdf tails and this property can only occur over 
limited interval of parameters. 

2 Chartists or Fundamentalists 

In paper I [T] and in [7] we have highlighted the key role played by the 
coexistence of two kind of agents (i.e. two kind of market strategies) 
and by the transitions from chartists to fundamentalists and viceversa. 
Anyway our model also shows some interesting aspects in the limit 
where the agent population is composed only by fundamentalists or 
only by chartists. For these limits we can derive some analitycal re- 
sults which are also useful to clarify various properties of the general 
case with both strategies. Morever this analysis will also suggest an 
approximated solution of the general case which will be discussed in 
section 5. 

Our discussion starts from the equation for price evolution (paper I [1] 
Eqs. (8) and (9)) that can be written as 

pt+i -pt = ^jfzj^Pt ~ p^^) + (1 ~ ^hipf - pt) + (^^t- (1) 

For the sake of simplicity we have dropped the explicit temporal de- 
pedence of the moving average pm and of x in Eq. ^ . 
In Eq. ([T|) the temporal scale is the elementary one which corresponds 
to a time increment At — 1. Our analitycal results will refer to this 
unitary time scale. On the other hand, computational results will be 
estimated from returns calculated each 100 time steps for convenience 
with respect to real data. Unless specifically stated the time scale of 
our discussion will the elementary one. The analitical formulas will be 
given in term of this elementary scale and in order to compare them 
with results on a scale which is k times the elementary increments, it 
will be sufficient to make the following substitution 

t ^ t' ^ kt. (2) 

In this way the new function of the new variable t' will be expressed 
in time units of length kAt. 

2.1 The case x = 

In this case Eq. ([T]) reduces to: 

Pt+l -Pt= lipf - Pt) + cr^t (3) 
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and, in order to avoid a trivial divergence oipt, we consider the interval 
< 7 < 1. 

Equation ([1]) describes the path of a random walker subject to a force 
which tends to revert her to pf and 7 is the parameter of the strength 
of this force. In fact the first term of the right member of Eq. ([T]) can 
be written as the derivative (changed in sign) of the quadratic potential 

In the following calcutation we set pf ~ wihout loss of generality 
since we are now considering the case of linear dynamics. Given the 
initial value po = (without loss of generality) , the explicit solution of 
Eq. ^ can be written at any time as 



Pt+i = aj2^,a-iy-'- (4) 

We obviously have tha10 E[pt\ ~ and the second moment of the 
random variable pt takes the form: 



It is sufficent to recall that 



(l-7)-2-l 



ifi^i 

1 if i = j 



(5) 



(6) 



in order to obtain Eq. 

Equation ^ shows that the second moment of pt is independent on 
time when i — > cx). In fact it can be shown that the auto-regressive 
process (AR) that defines the case x = is stationary. Hence in 
the following results of this section we assume that the limit to the 
stationary regime has already been performed so that the temporal 
dependence can be omitted. Furthermore it is worth noticing that 
time-dependent terms exponentially go to zero with speed fixed by 7 
which is typically 0.006 in our modeo 

We now turn our attention to the properties of the returns rt that we 
define as rt+A = Pt+A — Pt where A is the time lag at which we are 
studying the increments of the process (see appendix A for a discussion 
on the definition of the returns). From Eq. Q we find: 



n+A = cr(l - 7) 



t-1 



i+A-1 

E 



t-i 

j=0 



(7) 



and we can also compute the variance of returns (taking into account 
stationarity wa can write vt+A = fA)' 



^(1-7)- 



(l-(l-7)^) 



(l-7)-2-l 



(8) 



ets 



^with E[-] we indicate the expectation value of the expression inside the squared brack- 
^thus the process is already nearly stationary for t > 1000 
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that correctly reduces to the well known result for the RW (t^A when 
7 goes to zero. Equation ([8|) can be calculated in a similar way of Eq. 
(O by Eq. Q . The analysis of Eq. ^ with respect to A leads to 

f (J^ ff^^S', 7 A for small A 
I ^ %i^7i-i for large A 

We see that for small values of A the behavior of Eq. ([8]) is similar to 
one of the RW except for the size of fluctuations which arc reduced by 
the attractive term due to the pf. In the limit of large A, instead, the 
fluctuations tend to a constant differently from the RW. 
We now investigate the normalized autocovariancc function of returns 
Pr-A (■'") and of squared returns p^^ (r) as a measure of market efficiency 
and of volatility clustering respectively (they depend only on the time 
difference r because of the stationarity) . We recall that Pr^{T) and 
p^2^ (r) are defined as follows: 

PtI(t)^ — (10) 

PrAr) = (11) 

where cr^ = £^[7*^] and /i^ — E[rA]. The former one can be computed 
in a similar way of previous fomulas because all the terms arising from 
Eq. (Uni) can be reduced to the evaluation of expressions like one in Eq. 
(O. We find (it is a well-known results that the correlation function of 
an autoregressive model is characterized by an exponential decay): 

P.A(r)--^e^^""(-'')l^^<^-^e-. (12) 

Equation (fT2ll shows the case x = exhibits a negative correlation of 
price increments and this correlation has a slow decay rate since the 
characteristic time tq is about 7^^ « 167. This fact is not surprising 
because we have already noted that Eq. ^ describes a walker which 
is attracted by pf. To put it simply, if we observe an increment of a 
certain sign at time t, increments of the opposite sign are more likely 
observed in the following time steps beacuse of the constraint of pf. 
In a similar way we can find from Eq. (|lip : 

(r) = /(A, 1 - 7)e--l '"(i"^)! /(A, 1 - 7)6-^^^ (13) 

where /(A, 1 — 7) is a function of the parameters A and 1 — 7 whose 
explicit form is rather complex but not very relevant since it refers to 
the values of the prefactor of the correlation function. As expected 
the decay characteristic time of the autocovariancc function of squared 
returns is half the value found in Eq. (fT^ because the path of the price 
is a fractal with dimension 1/2. But, while a decay time of about 7"^ 
is rather long if we are studying the deviation from efficiency, a decay 
of 7~"'^/2 is very short with respect to the real volatility clustering. 
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Figure 1: The plot shows the comparison of Eq. [12] and Eq. [13] (solid lines) 
with the results of a simulation (o) and (□). We see that the autocovarinace 
function of returns, which is a proxy to test the market efficiency, is negative 
when X = 0, conversely the volatility has a positive correlation but it is 
negligeable. In fact both of autocovariance functions have a decay time 
which is or fraction of 7~^: in terms of efficiency this case is locally 
inefficient, instead in term of volatilty clustering this case substantially does 
not exhibit this phenomenon. 



From a financial point of view these two results mean that there is a 
predictability of price increments for about 200 price increments (at 
the elementary scale) corrisponding to a rather long inefficiency while 
the volatility clustering decays quicker that inefhency, so that the case 
a; = shows a behavior which is the opposite of reality. In Fig. [T]we 
compare our analytical results with computer simulations of the case 
a; = 0. Before proceeding to the next section it is interesting to make 
some comments on the continous limit of the case x = 0. In this limit 
the price equation becomes (p/ = 0): 

dp{t) = -jp{t)dt + adwt (14) 

where dwt is the differential of the Wiener process. The stochastic 
differential Eq. ([Tl)) is the well known Ornstein-Uhlenbeck process 
(for further details see |S]). The continous limit of Eq. ^ statisfying 
the same intial condition is: 

p{t)^(je-'^^ [ e'^dws (15) 

"'0 

where the integral of the right side of Eq. ([T5)) is a Ito stochastic in- 
tegral. In order to test the results of the discrete time case we can 
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compute the expressions corresponding to Eqs. ([8]) and (|12p in the 
continuous hmit. We have verified that the formulas of the discrete 
case correctly reduce to the continuous one in the appropiate limit of 
the parameters 



2.2 The case x = l 

The development of a bubble or a crash corresponds to a strong increase 
of the fraction of the chartists, here we are going to explore the limiting 
case in which all agents are chartists. In such a case .t = 1 Eq. ^ 
becomes: ^ 

pt+i ~ pt ^ j^—j{pt ~ Pm) + cr^t (16) 

This equation describes a random walk in a quadratic field potential 
as the previous case but this time the force is repulsive and the is- 
tantaneous (unstable) equilibrium position is the value of the moving 
average pm ■ 

This approach to model the stochastic process of the price was first 
introduced by [51 [TU] and following these papers we rewrite Eq. 
as a telescoping sum: 



Pt+i - Pt 



M 



(A/ + 1 - i){pt-i+i - pt-i) + cr^t 



(17) 



where a = b/{M{M — 1)). In this way we can find a simple recursive 
formula for a M-dimensional system whose first component is equiva- 
lent to Eq. ini): 

ft+i=Tft + Tft (18) 
where T is a Af x 71/ matrix defined as: 



/ 




02 . . 


■ ■ • ■ aM \ 
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. . 


. ... 







1 . . 


. ... 


V 





. . 


.10/ 



a{M+l-k) (19) 



and ft and fft are vectors M-dimcnsional so defined: 



rt = 



Pt -Pt-i 
Pt-i -Pt-2 



\ 



Vt = 



\ Pt-M+l - Pt-M / 







V / 



(20) 



The authors of [51 [TU] propose a simple approximation to calculate the 
variance of returns and find that the variance of the RW is magnified 
by a b-dependent factor larger than 1 so that the repulsive potential 
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magnifies fluctuations with respects to a RW. We note that it is possible 
to solve the vectorial Eq. (fT8|) recoursively: 



fe=0 fc=0 h=a k=l h=l 

where yo is the initial condition, for the sake of simplicity we assume 
that i/q ~ ioY i = 1, . . . , M . The solution of Eq. (fT6|) can easily be 
recovered by considering the first component of the vector rt- We 
find that E[r]] ~ and the variance is: 



fc=0 fc=l 



A-l 

(22) 



where xu^t - (E^To T*+''-'^)ii and pu = {Y.t=i T"^^)ii. The struc- 
ture of Eq. ([2H) is rather interesting, in fact it is composed of two 
sums: the first one depends explicitely on time t, while the second 
depends only on A. This is because it originates from the last term 
of the right member of Eq. ([2T|l that depends on t only through the 
index of the gaussian random variable Therefore this dependence 
will eventually introduces only a Kronecker's Delta function if some 
expectation values are calculated. In order to have a stationary state 
the Eq. ([^^ can no more depend on t and indeed in Fig. [2] we have 
verified that the process appears to converge to a stationary value for 
the variance when t grows. 

As counter check (and as a more formal test to verify the existence of 
the stationary state) we can evaluate numerically the spectrum of the 
matrix T. Naming Ai with i = 1, . . . , M the eigenvalues, we find that 
M(Ai) < 1 Vi and thus the time dependent term of Eq. (PT|) will behave 
like a geometric progression with ratio less than 1, that is the sums are 
convergent. So, as in case a; = 0, the process admits an asymptotic 
stationary state and in practice this means that the process is already 
stationary for t > 500. 

In Fig. [3] instead we compare our exact result with the approximation 
of the return variance given in [9l [1^ , the approximated formula gives 
a very good agreement in the asymptotic limit where the process is 
substantialy a RW. On the other hand the approximation does not 
reproduce the super diffusion for A < 1000. In fact at this time scale 
the repulsive role played by the moving average is crucial. 
In a similar way we can calculate the normalized autocovariance func- 
tion of returns Pr^ir): 

ELo(^M)^ + Et"i'(/3^-)^ 

The numerator of Eq. (j23p is positive-definite (and then p^A (''') too) 
because it is a linear combination of products of positive terms. In 
Fig. [¥]we report the estimation of Pr^ir) from a numerical simulation 
of the process in the case a; = 1 compared to the analitical result of 
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time 



Figure 2: In this picture we show the quick convergence to the stationary 
value of the variance of returns for time lags equal to 100. The existence 
of the stationary limit of Eq. ()22p is also supported by the fact that all 
eigenvalues of the matrix T are less than 1. 

Eq. (ESI). 

The autocovariancc of squared returns can be also computed in prin- 
ciple but its structure is very complex, in Fig. chartists the asterics 
corresponds to a simulation which shows its positive sign and the fast 
convergence to zero. Once again we have a process which is character- 
ized by a quite long inefficiency (positive this time), and on the other 
hand by the absence of persistency of volatility clustering. To sum- 
marize the studies of this section we see that the two extreme cases in 
which a type of agent dominates on the other one correponds to two 
different regimes. One of these are characterized either by gaussian 
fluctuations smaller than those observed in a pure RW and by nega- 
tive arbitrage (a; = 0). The other one instead has gaussian fluctuations 
larger than those observed in a pure RW and positive arbitrage {x = 1). 
Neither of these two regimes show volatility clustering, therefore the 
next question is how these ingredients are mixed when we take into 
account the dynamics of x so that it can assume values between and 
1. Next sections will answer this interesting question. 
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Figure 3: Eq. (j22p . which is the exact expression of the variance or returns 
in the case x = 1, predicts the superdifFusion of the process in the case 
X = 1 for small value of time lag A, in fact the slope of the red line is larger 
that the slope of the green solid line which behaves like a random walk. As 
observed in [9l [10] the approximation proposed by these authors works very 
well when A ^ 100 where the process is close to a random walk except for 
a magnification of the effective (7(0). 
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Figure 4: As in the case x = we show how and if the case x = 1 reproduces 
the Styhsed Facts, this case is characherized by a positive inefficiency whose 
temporal length is similar to the case of x = 0, instead once again we see that 
neither the moving average is able to reproduce a significant and persistent 
volatility clustering 
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3 DifFusional properties 



Diffusion processes are time-continuous Markov processes that have to 
satisfy three conditions (see [TT] for a mathematical formulation of this 
characterization). These conditions, roughly speaking, state that dif- 
fusion processes arc fully characterized by their drift (i.e. their mean) 
and by their diffusion coefficient (i.e. their variance). Diffusion pro- 
cesses are essentially gaussian processes which are fully specified by 
the knowledge of the mean and the varianc^. 

One should notice that these diffusional properties (limited to the 
squared returns) do not provide a complete information on the sys- 
tem in the presence of anomalous tails, however we are going to see 
that the diffusion analysis leads in any case to some interesting in- 
formation. Before moving to the statistical analysis we would like to 
point out that the process in the case a; = is gaussian and Marko- 
vian. Also the case a; = 1 is gaussian, being linear in the price and 
in the white noise, but no more Markovian due to the moving average 
which introduces a term that explicitly depends on the price path. The 
full process (Eq. ([T])) would be gaussian if x was fixed because of the 
linearity of Eq. ([T|) (but non Markovian). Nevertheless the random 
process which drives the switch of populations destroy the gaussianity 
and the return pdf actually has non gaussian tails. We will show this 
effect in section 5 where we propose an approximated mechanism to 
explain the non gaussian tails. 

It is well-known that the variance of the increments of a RW is pro- 
portional to the time scale A that is (j(A) = (TqA^''™ with ^„u = 0.5. 
We now investigate the scaling properties of the variance of increments 
of the process of Eq. H]) as a function of the time lag A at which returns 
are calculated for different values of N and we sometimes find devia- 
tions from ^rw We also investigate the variance of the increments on 
two different temporal windows, the former could be defined a short- 
scale analysis, the latter instead an analysis on long temporal horizons. 



Short time scale We report the plot of the variance pattern in 
Fig. [5] while in table [T] we show the corresponding results of a fit with 
a power law of the curves of Fig. [5] in the short scale caseQ- 
Once again we stress the crucial role played by the number of agents N 
with respect to the statistical properties of the model. In fact we clearly 
sec in table [T] that the variance, for a fixed A, is a decreasing function 
of A^. This can be understood in term of the population dynamics (see 
section 1 and [TJ[7]: when the number of agents N increases, the frac- 
tion X of chartists tends to zero because of the exponential term of the 

^however a gaussian process can be non-Markov differently from a diffusion process 
which is necessary a Markov process 

^tlie evaluation of the error bars when we deal with dependent variables is longly 
debated so we prefer dropping them and consequently we propose the fitted parameters 
without indication of their statistical errors 
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asymptopic stable distribution of Eq. (21) of paper I [T] and chartists 
have previously been recognized as fluctuation magnifiers, while fun- 
damentalists as fluctuation reducers. 

The most striking result is however the value of jU found for N « 500, 
when it is nearly equal to ^rw ■ The increments have the same scaling 
behavior of a RW but it would be misleading to conclude our pro- 
cess is a RW because we know that the process is non gaussian and 
non Markov since it exhibits non trivial correlations and non gaussian 
tails. This implies that a complex intermittent behavior may not be 
detected by a simple diffusion analysis. It is worth noting that the be- 
havior for any N is bounded by the limitating cases x = and x = 1. 
The apparent RW-like scaling behavior which for our parameters occur 
around N = 500, is the results of the a balance of the opposite ten- 
dencies due to the fundamentalists and chartists. A similar balancing 
situation has been discussed in [TU] even though the framework is 
rather different than the present one. 
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Figure 5: Investigation of the scaling behavior of the second moment of the 
returns, that is the diffusional properties, for x = 0; 1 and for different values 
of N and in the range of small values of A. In the main plot we see that the 
short scale behavior is approximately of the type and the exponent fi 
varies from the maximum 0.7 when x = 1 (superdiffusion) to the minimum 
0.44 when x = (subdiffusion). When N ~ 500 the behavior converges to 
one of a random walk. The small insert shows instead the begining of the 
regime dominated by fundamentalists which develops in the limit of large 
A. 



N 




N 




RW 


0.500 


300 


0.523 


x = 


0.439 


400 


0.509 


X = 1 


0.701 


500 


0.496 


10 


0.609 


600 


0.482 


50 


0.596 


750 


0.466 


100 


0.581 


1000 


0.455 


200 


0.555 


5000 


0.440 


250 


0.541 







Table 1: Results of fit of the exponent for A £ [1, 100] 
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Long time scale When wc turn our attention to a longer temporal 
horizon (A ^ 100), wc can observe very different features (see Fig. [6]). 
According to Eq. ^ the variance of the case a; = goes to a constant 
when A — > oo while from Eq. ^ the variance of a: = 1 grows as 
in a RW when A ^ oo {fJ-x=i is larger than 0.5 at short times and 
asymptotically tends to 0.5 from above). 

The balancing effect fails in this limit (in Fig. [6] we show only the case 
N = 500 because for all N the variance tends to a constant but this 
constant value depends on TV). In order to explain what happens when 
A grows we observe that the characteristic time scale of the moving 
average is M (or few times M). Therefore, on time-scale larger than M, 
a coarse-grained version of the process relate to would be nearly 
undistinguishable from a RW. Instead the case a; = does not exhibit 
diffusional properties since 

Coo = lim = (24) 

where Coo is the asymptotic diffusion coefficient. In this way if A 
becomes very large (i.e. the coarse-grained limit) the term relative to 
the moving average behaves like an additional noise. In fact, as we 
have seen in Fig. wc just observe a greater effective variance than 
the one predicted by Eq. fH]). 
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Figure 6: Scaling behavior of the second moment of returns. Diffusional 
properties for x = 0; 1 and for the full mixed process with = 500 in the 
limit A — > oo. While at the short time scale the subdiffusion of the case 
X = and the superdifFusion of the case x = 1 someway balance out, when 
A grows the variance tends to a constant as in the case x = 0. 



4 Transition to gaussianity? 

In the previous section we have studied the A— dependence of the vari- 
ance, but the parameter A may a priori play a role in the sliape of the 
return pdf too. The motivation of such an interest can be traced in 
many analysis of market data (see [TH [T3]) where data appear to be 
well- fitted by a gaussian when A oo. This phenomenon is still un- 
clear, some authors invoke the central limit theorem (CLT) so that the 
validity of CLT, that fails for small value of A, seems to be recovered in 
this limit although it is still unclear if CLT or critical-phenomena-like 
arguments should be used or neither of them. 

We are going to discuss in detail the problem of the Fat Tails within 
our model. In Fig. [7] we display the pdf of the returns El while in Fig. 
[5] we show the same results but we rescale the returns dividing them 
by their variance. From these figures it is evident that in our model 
Fat Tails do not dissapear for A ^ oo, this may appear contradictory 
with the results of Fig. [5] where the variance is shown to converge to 
a finite value for A — > cx). We arc going to see that the origin of this 
apparent discrepancy is due to the fact that the second moment is not 
appropriate to discuss the Fat Tails because it is dominated by the 

^data are from a very long simulation t = 10® 
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central part of the pdf. 

In Fig. [7]and allpd2 we see that all the pdfs have the same behavior: in 
the central region they appear to follow the shape of a gaussian while, 
in the region of rare events, pdfs have tails decaying slower than gaus- 
sian. Therefore the interaction between chartists and fundamentalists 
must be the key element to increase the probability of observing large 
fluctuations. An analysis of the data shows that the central gaussian 
behavior is essential due to the action of fundamentalists {x ~ 0) as 
discussed before. 

The fact that the rescaled return pdfs do not collapse everywhere in 
the same curve confirms that the scaling property of the variance does 
not fully characterize the process and thus the rescaled returns exhibits 
only a partial collapse that indeed fails in the tail region (see Appendix 
C). 
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Figure 7: Probability density function of returns for different values of time 
lag A at which returns are performed. Our minimal model does not seem to 
show a transition to a gaussian regime since the non gaussian tails persists 
even for very large values of A. 



In order to answer why in this mininal ABM deviations from gaussian- 
ity do not disappear when A — s- cx) we have to consider the microscopic 
mechanism which leads to the tails. 

On this account Fig. [9] is very instructive: we show the return pdf 
for N = 500 and the corrisponding pdf of the cases x = 0; 1 for the 
same A. We see that the central part of the curve for N = 500 (blue) 
is very close to the gaussian of the case a; = 0, while the Fat Tails 
appear to be somehow caused by the increasing width of the gaussian 
of the case x ~ 1. Since we have already seen that the model would 
be gaussian if x was fixed, we can conclude the non-linearity deriving 
from the dynamics of x must be the key to understand the emergence 
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Figure 8: Probability density function of returns normalized by their vari- 
ance. The non collapsing tails suggest that the scaling behavior of the 
variance captures only partially the process. We will try a different normal- 
ization of returns in Appendix C but we will not obtain collapsing tails for 
any A. These difficulties arise from the fact that tails draw origin from the 
non linearity introduced by the population dynamics. 



of non gaussian tails. 

The possible relation of our results to real data should be considered 
with great care. For example our model is tuned to elucidate the ori- 
gin of the SF and for simplicity we have set pf to be a fixed value 
independent on time. This is useful from a conceptual point of view 
but it is certainly non realistic assumption. As we arc going to see in 
section 6 the introduction of a RW dynamics for pf{t) would change 
completely the asymptotic behavior of the variance and would lead to 
a real crossover to gaussian behavior at long time. 



5 Nature of Fat Tails 

We now propose a simple approximation which clarifies the origin of 
the non gaussian tails. In section 2 we have found that the return pdf 
of the case x = and x ~ 1 are respectively 

2 

gf{y):=g(y\x^O)^^L=e^^ (25) 



gM := g{y\x = 1) = —===e (26) 
v/27rcr2(A) 

with a J Ri cost when A ^ I and crc(^) ~ "'co^- other words 
we know the solution of our model in this two limiting cases. It is 
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Figure 9: The four panels show the comparison of the return probabiUty 
density function for the extreme cases a; = 0; 1 (which are gaussians) and 
for the mixed dynamics for = 500 for different values of A. The growing 
gaussian of the case x = 1 seems to stretch the non gaussian tails which 
are produced by a non trivial interaction of these two gaussians through the 
population dynamics, that is through by the temporal dependence of x. 
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instructive to make the simple assumption that the general solution 
r^v for a certain N can be written as a superposition of the two cases 



rN = xTc + (1 - x)rf (27) 

where rf,rc are the returns of the case x ~ and x ^ 1 distributed 
respectively as Eq. and Eq. PH|) and we assume stationarity. 

Hence the pdf of r^, for a fixed valaue of x, is given by the following 
convolution 

p{r\x) ^ J J dydzgcjiy, z) A{r - xy - (1 - x)z) (28) 

where gcj{y, z) is the joint pdf of chartists and fundamentalists. 
We neglect the interactions between chartists and fundamentalists and 
in this approximation we can factorize gcj{y,z) = g f{y)gc{z) and in 
the end we find that 

P{r\x) = -=L=e-^ (29) 
V 27rcr^ 

where = x^a^{A) + (1 — x)'^a'j. This result is trivial because the 
convolution of two gaussian is a gaussian. 

However, a much more interesting result can be obtained by consider- 
ing the equilibrium distribution feq{x) of cc, as given by Eq. (21) of I 
[T] and integrated Eq. ((^ over this distribution of values of x. This 
integral can be numerically evaluated 

E,[p{r\x)]^ f dxfeg{x)p{r\x). (30) 
Jo 

In Fig. \W\ we compare the approximated curve corresponding to Eq. 
(PU]) with the simulated one and we can see that there is a good agree- 
ment especially for A < 10'''. 

From Fig. [TO] one can see that, at short time, the superposition 
(A < 1000) is a very good approximation, but for longer time this 
approximation leads to an over-estimation of the tails. This can be 
undestood rather easily by considering, for example, the behavior dis- 
cussed in relation of Fig. [6l For long times the variance of the chartists 
alone diverges as in a normal diffusion process, on the other hand the 
variance of the fondamentalists alone is bounded to a fixed value. Con- 
sidering the attractive role which pf plays in dynamics of the complete 
model it is clear that the variance of the complete model will not di- 
verge. In this respect the incoherent superposition of the two limiting 
cases leads to a excess of variance with respect to the real dynamics of 
the model. 



This simple superposition approximation can also be used to test which 
features of the model depend crucially on the chartist-fundamentalist 
interaction. Assuming that the general solution can be put in the form 
of Eq. (|27p and assuming that Tc and rj are indenpendent and sta- 
tionary processes, we can compute the approximated autocovariance 
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Figure 10: Comparison of the approximation proposed in Eq. ()30p with 
the probabihty density function of returns estimated from a simulation for 
different values of A. The approximation gives a good agreement with sim- 
ulation data until A < 1000. For larger values of time lag A the tails are 
overestimated by the approximation. This corresponds to an overestimation 
of the effective fraction of chartists when the average on x is performed. This 
is in agreement with the behavior of the variance shown in Fig. [6l 
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function of returns (r) and of squared returns p^2 (r) . The former 



IS 



2 ^2 



Pr-i {r) + 



Prii^) 

(31) 



where cr? and p / (r) are respectively the variance and the normahzed 

autocovariance of returns for x = and cr^ and pr<= (''') are the same 

quantities for x = I. 

In the same way the latter is 



Prl{T\x) = 



{l^x^ Mj (l-x-) aio] 
x^ M4 a;2 M5 



6- 



t2^2 



(l-a;)4M* {\-xY Mf 
x^ M4 (1 - x)2 Af4 n -1 



6(1 — x)2 0-2^^2 g2;2 ^2^2 



P{rl)^iT) + 

-1 

Piriy-ir)+ 
Pr^Jr)p^fiT) (32) 



where Alj and are the fourth moment of returns in the case x = 
and a; = 1. 

As in the previous section we have to evaluate the following expectation 
values to elimate the dependence on x in Eq. ([31]) and Eq. ([32]) 



(33) 
(34) 



Pr^i^) = I dx Pr^{T\x)Jeq(,x) 



PrliT) ^ dxPr2 {T\x)feqix). 

Jo 



Fig. [TT] shows the comparison of these two autocovariance functions 
calculated from a simulation with the corresponding formulas just 
found. The superposition approximation works quite well for the au- 
tocorrelation of the returns but it fails to reproduce the correlations 
of the squared returns. For the simple returns we conjecture that this 
agreement is due to the linear balancing of the positive inefficiency 
of chartists and the negative one of fundamentalists and corresponds 
really to the incoherent superposition of opposite effects. For volatil- 
ity clustering instead the simple superposition dose not capture at all 
the nature of the phenomemon. The coherent persistence of a chartist 
bubble, in fact, cannot be reproduced by the simple superposition ap- 
proximation. On the other hand volatility clustering are not present 
in any of the two limitating cases. 
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Figure 11: The simple approximation proposed in this section also allows to 
extimate the autocovariance of returns and of squared returns. The global 
efficiency of the model seems to be well reproduced by the simple super- 
position of the solution of the case x = 0; 1. The approximation fails in 
reproducing volatility clustering (even for t < 1000) because the individual 
cases X = 0; 1 do not reproduce this Stylised Fact. This comparison clari- 
fies the key role of the population dynamics (i.e. of x) in order to generate 
bubbles which gives non-zero correlation for volatility. 



23 



6 Random walk of fundamental price and 
effective gaussian transition 

We have seen that our minimal model, strictly speaking, does not 
lead to gaussian transition for the returns at long times. The fact 
that real data seem to show such a transition therefore may appear as 
problematic. On the other hand the model was aimed at the identi- 
fication of the nature of the SF and we made the simple assumption 
that Pf = cost. Clearly this is unrealistic and for a real comparison 
with the data one should consider also the dynamics of Pf{t). This 
will change the situation in an important way and will reconcile this 
apparent discrepancy. 

Following this reasoning we could expect that the return pdf gets very 
similar to the pdf of the process of Pf{t) when the term due to the 
drift becomes leading. Hence a natural choice for the pf consists in a 
random walk with a variance smaller than the variance of the white 
noise of Eq. ([T]). How much it is smaller just fixes the time scale A* 
at which the drift becomes the leading term. 

In Fig. [12] we report the same diffusional analysis of Fig. [6] in the case 
of (Jp^ ~ lO^^CT where a is the same parameter of Eq. [T]and Fig. [T3] 
shows the pdf for various value of A. The main difference with Fig. 
[7] is represented by the brown square curve which is well-fitted by a 
gaussian without appreciable tails. This results is obvious because we 
have choosen that pf makes a random walk and the brown curve is 
in the regime A > A*. The random walk of pf makes the variance a 
strictly increasing function of A (for A > A*) and the data in Fig. [T?] 
provide a clear evidence that the dissappearence of the non gaussian 
tails is due to a statistically irrelevance of them when an additional 
gaussian term is present. In this respect one cannot really say that 
Fat Tails disappear but they become extremely rare events and are 
essentially obscured by the gaussian behavior of pf. 
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Figure 12: We show the same analysis of Fig. [6] for the scahng behavior 
of the variance of the returns in the case in which the fundamental price 
is not constant (as in Fig. diffulong) but it follows an intrinsic random 
walk with apf = 1/10). We note, as expected, that when A — > cx) the 
leading contribution comes from the random walk of p/. Therefore in our 
interpretation the transition to gaussianity of real data at long times arises 
from the predominance of the pj random walk with respect to the Stylesd 
Facts which become extremely rare events at long times. 
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Figure 13: The panel shows the comparison of the return probability density 
function for different A. When A = 10^ this function is well fitted by 
a gaussian (solid line). The effective transition to the gaussian behavior 
derives from the random walk oipf{t). 
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Figure 14: The random walk of pf makes the return probabihty density 
function of the case x = 1 wider and wider therefore, when N = 500, the 
deviation from the density of the case a; = is not observed. 
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7 Multiplicative version 



The role of Walras' law (see [14]) in economics is very similar to the 
role of Newton's law of the dynamics in physics because Walras' law 
introduces a relationship between the excess demand (ED) and the con- 
sequent variation of the price that ED causes (the excess demand is 
usually defined as a function of the unbalance of demand and supply) . 
This law has been formalized in many ways since it was introduced 
about two centuries ago and the version, which belongs to the frame- 
work of our model, is the following 

-^^PED (35) 
p at 

or when time is discrete (At = 1): 

^^[3 ED (36) 
P 

The basic idea of Eq. (|35p is that price follows a multiplicative process 
rather than a linear one so its percentage increments, rather than the 
increments themselves, are proportional to ED. There is appreciable 
evidences in the experimental data that a multiplicative dynamics is 
indeed more appropriate than a linear onc [15l ITC] , 
Up to now we have neglected this multiplicative nature in our model 
and we have adopted the simplified linear dynamics. As we are going 
to see in this section the two cases essentially coincide for small fluctu- 
ations. The purpose of our model was to trace back the detailed origin 
of the SF and their self-organization. In this perspective the linear dy- 
namics does not change the essence of these elements but it represents 
an useful simplication especially in deriving some analitical results. In 
this section we discuss in detail the analogies and differences of the two 
cases. 

In our model the (linear) price increments are proportional to the term 
~xb/{M—l){pt —phi) + {l—x)"f{pj—pt) which we identify as the ED, 
in this respect it corresponds to a linearized Walras' Law. Considering 
this as the small price increments limit of a multiplicative dynamics wc 
can now try to go back and consider the real dynamics from which this 
approximation may arise. Wc arc going to sec that this backward re- 
construction from linear to multiplicative presents some subtle points 
which are worth to discuss in some details □ 

7.1 Choice of excess demand 

Wc would like now to construct a price equation which is consistent 
with Eq. (|36p so we define the price increments as a geometric random 
walk with an additional term that is indeed the excess demand: 

ln(pt+i) = In(pt) + ED + a^t (37) 



in section 7 and in all its subsections returns are defined as rt+A = \n{pt+A/p. 
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where are independent and normally distributed random variables. 
However Walras' law does not specify which is the correct choice for 
ED, apart from the fact that the dimensional analysis of Eq. ([55)) 
implies that ED must be somehow normalized (dimensionless). The 
term of Eq. ([T]) that plays the role of ED has the dimension of a price 
so the simpliest way to make the ED of the linear version dimensionless 
is to normalize it with the quantity which has the dimension of a price. 
This normalization is not univocous and for example one could have 
the two different situations 

ED = (^^LZM + (1 _ (38) 

M -I pt pt 

ED = ^^^ll^lEl + (1 _ (39) 
M - 1 pm Pf 

with obviously pf ^ 0. 

In Eq. psp chartists and fundamentalists have the same benchmark 
to calculate the percentage ED while in Eq. ([55]) each of them has 
different one: chartists use the moving average, fundamentalists use 
the fundamental price. In Fig. [15] we show and compare the results in 
the two cases. The behavior is very similar and the detectable feature 
is an increase of the fluctuations in the case of Eq. ([5^ . 



Despite these similarities the two normalizations correspond to a quali- 
tatively different nature of the dynamics. In fact suppose, for example, 
that the price has just reached the local minimum. Then it begins to 
grow and so it crosses the moving average from below. Once crossed, 
the following inequality holds (p4—pj\/)/PM > {Pt^PAi)/pt- Conversely 
if the price has just reached the local maximum the opposite inequality 
is {pt ~ pm)/pm < {pt ^ PM)/Pt- This situation would imply that in 
the limit dominated by the chartists the price will unavoidably go to 
infinity for Eq. ([55]) or to zero for Eq. ([55]) depending on the initial 
condition. This simple example shows that the two normalizations can 
affect in a very different way the large fluctuations of the dynamics. In 
our opinion the most appropriate normalization is the one correspond- 
ing to Eq. ([38p because pt is a natural value recognized by all agents 
and in addition this leads to the correct linearized dynamics as now 
we are going to see. 



7.2 Calibration of parameters and model instabili- 
ties 

Most of the ABMs share the common aspect that they (more or less) 
reproduce the empirical facts (i.e. SF) only in a certain region of the 
phase space of their parameters. A tipical example is the model of Lux 
and Marchesi [2 [3] or [U [7] . Therefore they require a fine tuning of all 
the parameters and even of the total number of agents N in order to 
reproduce the SF. 

This situation is already present in the models with the linear dy- 
namics. The multiplicative version, as for example in our model Eq. 
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Figure 15: The two different normalization for ED lead to similar Stylised 
Facts. In our opinion (see text) the most natural choice is to normaHze ED 
with the price itself which is an objective piece of information shared by all 
agents. 
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(|37p . is dramatically unstable and more sensitive to small changes of 
parameters than the linear one. The results is that the dynamics of 
price cannot be simulated for a large part of the phase space of the 
parameters because it is subject to extreme singularities towards zero 
or infinity. Consequently the main attempt of tuning the parameters 
consists in finding the subtle equilibrium between the weight of the 
noise (its variance a) and the weight of the ED (6 and 7) in the evolu- 
tive Eq. ([571) of the price. 

The starting point of this non trivial calibration is the choice of a 
reasonable value of the variance of the white noise. In fact if the vari- 
ance a would be of order 1 the typical increment in a single time step 
would be about 100% of the price or more, which is totally unreason- 
able. On the other hand, in order to calibrate the noise we should 
relate the model time step to real time. The time scale of our model is 
expressed in arbitrary units and these should be tuned with respect to 
the realistic dynamics. If we would like to intepret our model time step 
as a single elementary operation in the real market this can be done 
in the following way. The typical price fluctuations of the price of a 
stock during one day is of the order of 1-2 %. The number of individual 
operations for a typical stock is of the order of 1000 — 5000. From these 
heuristic arguments one can the argue that a typical variance in the 
model should be of the order a = 0.001. 

Then we should set the parameters b and 7 so that the magnitude of 
ED is comparable with the white noise. In fact if it would happen 
< ED^ the process will be dominated by the gaussian noise 

and it would not show the non gaussian tails. In the other hand if the 
weight of ED is too large we will observe the intrinsic unstability of the 
process due to the exponential dependence on ED in Eq. (|37|) . Wc find 
that the upper values that can be choosen are b w 1.7 and 7 w 0.01. In 
Fig. [iniwe compare the tails of the return pdf for two set of parameters 
{b = 1.7, 7 = 0.01, a ^ 0.00112} and {6 = 1.0, 7 = 0.006, a = 0.00112}. 
It is important to remark that the former set gives a significative devi- 
ation from the gaussianity in the sense of the SF, while the latter set 
instead produces a pdf which is very close the original gaussian. 

We have argued that the linear version should be a linearization of 
the original Walras' law Eq. ([55]) . If we now linearize Eq. ([57]) with 
the normalization of Eq. (jM]) , it follows that 

Pt+i -~pt~ {^Pt)^t + xjj^—^{pt -Pm) + (1 - a;)7(p/ ~pt). (40) 

This linearization is justified if (pt+i — Pt)/pt 1 and this condition 
is generally true for our choice of parameters. The comparison of Eq. 
(fin)) with Eq. ((T)) highlights that the only difference lies in the variance 
of the white noise. This variance is a random variable in Eq. (|40p and 
the linear version Eq. ([1]) is essentially the linearization of the multi- 
plicative dynamics if the variance of the white noise can be assumed 
as nearly constant and its stochastic nature negletcted. Let us define 
the effective variance as tr = apt, the random variables and pt are 
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Figure 16: The multiplicative process results appear to be extremely sensi- 
tive of the region of parameters choosen. Here we show two slightly different 
set of parameters which lead to drastically different results with respect to 
the Fat Tails phenomenology. 

independent beacuse pt = Pti£,o,S,i, ■ ■ ■ jCt-i)j hence we can focus just 
on the effective variance a. In the average E[a] = Pfcr and the typical 
fluctuation f2, normalized with the squared average, is 

The linear process of Eq. ^ is therefore justified if the fluctuations 
given by Eq. (j¥T|) are negligible. In other words Eq. ([T]) is the limit 
of the linearized Eq. (|40|) when goes to zero. In table. [2]wc show 
that fluctuations around the mean value of a are very small indepen- 
dently on pf. This shows that for a choice of the parameters which 
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Table 2: Estimation of for different value of pf. 
appears reasonable to reproduce real markets this effect is only of the 
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order of 3%. Therefore the hnearized version of the dynamics indeed 
is a vahd approximation to the multiphcative one. It is important to 
note however that a different choice of the parameters could lead to 
a situation in which even in the limit of small price increments the 
simple linearized (with the random noise) form would be incorrect. 
Unfortunately we are not able to perform an analitical approach for 
the multiplicative version as we have done in section 2 because of two 
reasons: the first one is the non linearity of Eq. (|57|) which makes the 
problem analitically very difficult, the second one is that the process 
cannot be simulated for x = 1 because the price goes quickly to zero 
or grows indefinitely according to the normalization used. 

7.3 Stylized Facts 

The previous discussion may suggest that the multiplicative version 
introduces only disadvantages because it is much more unstable in 
comparison with the linear version with respect of phase space of the 
parameters. The limit x = 1 does not exist and even the limit x = 
appears impossible to be studied analitically. In general very few 
results can be derived analitically for the complete model. But we 
have also important advantages, the first one is that the multiplicative 
dynamics is more realistic. We are going to see in fact that the SF of the 
multiplicative version are substantially more pronounced with respect 
to those of the linear version. In particular the multiplicative dynamics 
shows an interesting relation between the power of returns and their 
degree of correlation, which is in better agreement with real data. The 
specific results we are going to show corresponds to the normalization 
of Eq. (PSj) . Even though stationarity is certainly weaker than in the 
linear dynamics also the multiplicative case appears to converge to a 
quasi-stationary state but with larger fluctuations. 
In Fig. [TT] we show the return pdf for the linear version and for the 
multiplicative one. The non gaussian tails are much more pronounced 
in the multiplicative case. In order to compare the two cases we have 
rescaled the returns normalizing with their variance. 

In Fig. [iH] we compare the normalized autocovariance of returns and 
squared returns for the two cases. The two models give similar results 
from a qualitative point of view but the degree of correlation of squared 
returns is nearly twice in the multiplicative case. Volatility clustering 
is therefore enhanced by the multiplicative dynamics. 

We can also investigate the degree of correlation of |r|'^ as a function 
of (f). It is well-known that the largest correlation in empirical data is 
observed for = 1 with respect to higher powers ([H]). In Fig. ?? 
and[2n]we report the autocovariance functions for different values of 4>- 
The linear version exhibits a maximum in correlation for = 2 while 
for the real data (see insert of Fig. ??) the correlation is larger for 
(/) = !. In Fig. [20] we can see that the multiplicative version leads to 
a realistic order of the degree of correlation with the maximum 0=1 
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Figure 17: Fat Tails for linear and multiplicative dynamics (time interval 
equals 100 steps). The multiplicative case is characterized by much stronger 
non gaussian tails (the returns are normalized with their variance) . The solid 
line represents the gaussian benchmark. The small panel on the right shows 
the non normalized returns probability density function in the multiplicative 
case. 



(consistent with market data). 

At the moment we do not have a full understanding of the origin of 
these different behaviors for the two dynamics. A naif intepretation 
could be to argue that in a multiplicative model returns are generally 
smaller than 1, differently from the returns in the linear version. On 
the other hand that our autocovariance functions are normalized so 
that they are correctly rescaled to be compared. A simple simulation 
of the linear version with a set of paramaters so that \rt \ < 1 for most 
of the time confirms that this arguments is doubtful because the same 
anomalous order in correlation functions is observed. 
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Figure 18: Volatility clustering and market efficiency for linear and multi- 
plicative dynamics. The autocovariance function of returns and of squared 
returns are qualitatively similar in the two cases. The volatility clustering is 
larger in the multiplicative case than in the linear one. The linear correlation 
shows very similar behavior. Note that in this minimal version of the ABM 
both the linear and multiplicative dynamics lead to a single characteristic 
time scale and therefore volatility clustering has an exponential decay. The 
possibilities of multiple time scales for the agents can modify this behavior 
in a more realistic one (power law like) as shown in Fig. 7 of [1] 
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Time lag (100) 

Figure 19: Linear model: degree of correlation of absolute returns with 
respect their power (j). The linear dynamics leads to a maximum correlation 
for (j) = 2 while the real data (see insert) show a maximum correlation for 
(j) = 1. The real data corresponds to daily dynamics of typical NYSE stocks 
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Figure 20: Multiplicative dynamics. In this case the order of strength of 
correlation functions is in agreement with real data corresponding to the 
maximum for (j) = 1. 



7.4 Diffusional properties and transition to gaus- 
sianity 

The second order properties of the multiplicative process and the scal- 
ing behavior of the variance of return pdf can be studied in the same 
way as in section 3 and 4 and all the considerations about the statisti- 
cal meaning of this analysis are the same except for the fact that the 
multiplicative process docs not allow us to simulate the case x ^ I 
(only chartists). 

The behavior of the variance of the increments for small values of the 
time lag A is plotted in Fig. (PT|) while tabic [3] is the analogue of 
table [1] where we report the results of a fit of the scaling exponent for 
different N. 

We obtain the same results of section 3 in which the process has a 
variance which scales appoximately as a RW for almost any value of 
N. In the same spirit of Fig. E] we study the scaling behavior of the 
variance on a very long temporal horizon and we find that for every N 
the variance goes to a constant (which depends on TV) when A — *■ oo. 
The value of this constant values is always larger than the asymptotic 
value of the case x = 0. 

In a similar way we can also investigate the dependence of the shape 
of the returns pdf on the time lag A. This analysis is reported in Figs. 
[22l and23l Once again the multiplicative version confirms our previous 
observation of persistent non gaussian tails for all values of A. 
The general situation is therefore extremely similar for the multiplica- 
tive and linear case. Consequently we expect that also in the multi- 
plicative case the introduction of a RW dynamics for p f (t) would lead 
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Figure 21: The variance of the returns for the multiphcative case is similar to 
the one observed in the hnear case. The diffusion exponent is approximately 
a decreasing function of A^. 



to an effective transition to the gaussian behavior without Fat Tails 
analagous to the one discussed in section 6. 



N 




N 




RW 


0.500 


300 


0.552 


X = 


0.487 


400 


0.535 


X = 1 


/ 


500 


0.523 


10 


0.632 


1000 


0.492 


100 


0.606 


5000 


0.489 


250 


0.561 







Table 3: Results of the exponent fj, for A € [1, 10] 
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Figure 22: Return probability density function for different values of A 
in the multiplicative case. Once again the tails do not disappear and the 
transition to gaussianity is not observed. Unless one introduces a gaussian 
dynamics for Pf{t). 




Figure 23: The normalized return probability density function shows a bet- 
ter collapse in a universal curve in the multiplicative case with respect to 
the linear one. However, since in this case the case x = 1 does not exist we 
are not able to make a discussion similar to one of Fig. El 
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7.5 Self-Organized Intermittency in the multiplica- 
tive dynamics 

Considering that multiplicative case shares the same population dy- 
namics and a very similar behavior with respect to with the linear 
version, we can consider the problem of the self-organization along the 
same lines. Thus the key point is again to allow that N can be, in its 
turn, a variable which is not necessarily fixed a priori. The basic idea 
is that an agent generally perceives small fluctuations to be not very 
interesting to trade while large fluctuations can be caught as the sign 
of trading opportunities. As in [7] we introduce an indicator a{t,T) of 
fluctuations which agents consider. 

<j{t,T) = ^ E (P^-P)' (42) 

i=t-T 

The agents enter or exit from the market depending on whether these 
fluctuations are larger or smaller than some thresholds 

a(t,T)>e„, (43) 
a{t,T)<Qout- (44) 

Clearly the thresholds 8i„ and Qout are different from those in the 
linear case and should be suitably adjusted to the properties of the 
multiplicative dynamics. 

We report in Figs. [M] the same analysis for the multiplicative case 
of Fig. 4 of [3 for the linear case. In full analogy with the linear 
case we can see that also in the multiplicative dynamics the volatility 
fluctuations decrease by increasing and only in an intermediate value 
of N {N w 500) corresponds to the intermittent dynamics which leads 
to the SF. This situation, toghethet with the threshold of Eq. |43] and 
Eq. 131] leads to the self-organization of the system around the quasi- 
critical value of N (N*), independently on the initual values N{t = 0). 
In fact if we start from a small value of N (i.e. N ^ N*) the large 
fluctuations will attract new incoming agents. Conversely if ^ N* 
the typical small fluctuations makes the number of agents decrease. 
When N reaches N* the system is relatively stable: the value of N* 
depends on the the choice of the parameters, mainly on the thresholds 
and it corrcspons to the intermittent state with the SF. 
The overall picture is very similar to the linear case with the only 
difference that the convergency to the quasi-critical is somewhat slower 
with respect to the linear case for large N. 
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time 

Figure 24: Volatility fluctuations for various vlaues of A^. As in pQITj we see 
that for large values of N (N = 5000), the volatility is very small and this 
situation is not interesting for the agents. For small values of N (N=50), 
instead, the volatility is always very high and the market offers arbitrage 
opportunities for the agents. If = 500 we observe an intermittent behavior 
as in real market. The insert reports the histogram of the main plot. 
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Figure 25: Self-Organized Intermittency. The mechanism of self- 
organization is the same of the linear case proposed in [H [7]. When N is 
larger than N* (N* ~ 500 is the number of agents that shows SF) fluctua- 
tions are typically lower than 0j„ therefore agents leave the market. Instead 
when N is smaller than A^* the fluctuations are larger than Qout and the 
agents enter into the market. 
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8 Summary and perspectives 



In this second paper we have considered many statistical properties of 
the minimal ABM introduced in [1] [7] . We have also considered the 
generalization of the model to the more realistic case of multiplicative 
dynamics. 

The model was introduced to focus on the minimal elements which are 
able to reproduce the main SF observed in real markets and to discuss 
their self-organization. These elements are essentially: chartist agents, 
fundamentalists agents, herding dynamics and price behavior. From 
the studies we present, we can make some comments about the role of 
each of these elements. For example, the fundamentalist-chartist com- 
petition corresponds essentially to the competition between stabilizing 
and destabilizing tendencies on the price. We believe that our specific 
assumptions to describe these tendencies are not crucial for the general 
behavior and that the results of the model would be robust as soon 
as these destabilizing and stabilizing tendencies arc present in some 
form. On the other hand the herding effect is absolutely crucial and 
represents the key element for the dynamics of the Fat Tails and their 
self-organization in this class of models. Therefore it would be very 
important to find a way to test the concept of herding in real markets. 
The specific price behavior in terms of F and C strategies is important 
in defining the general properties of the strategies (stability versus in- 
stability) but it is not crucial with respect to the specific nature of the 
SF. 

We have mostly considered the model within its linear dynamics be- 
cause it represents a powerful simplification in view of analitical re- 
sults and interpretation of the data. The more general multiplicative 
dynamics is shown to be extremely more sensitive to the parameter 
range and much more difficult to study in all respects. However, the 
main results for the origin of the SF and the their self-organization, 
which we have obtained from the linear dynamics, have been also con- 
firmed by the non linear dynamics. The main difference is that the 
non linear dynamics enhances the deviation form the gaussian behav- 
ior and makes the SF more marked. 

The phenomenon of self-organization intcrpretated as Self Organized 
Intermittency (SOI) and Finite Size Effects are fully confirmed also by 
the non linear dynamics. 

Our minimal ABM can be easily generalized to consider realistic el- 
ements which are important in the real data but that are not con- 
sidered crucial (in our perspective) to clarify the origin and the self- 
organization of the SF. One of these elements is the value of the fun- 
damental price p f which we consider as constant (or zero in the linear 
model). It is easy to introduce a suitable RW for Pf{t) which permits 
to understand the dissapearance of the SF and the transition to the 
gaussian behavior which is observed in real data at long times. 
It is easy to consider other possible improvements to make the model 
more realistic. Given the microscopic understanding we have achieved 
these can be now analysed one by one in a systematic way. For exam- 
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pic an clement which is considered important [17|, I18j and which we 
intend to study in the future is the problem of the finite liquidity in 
the market. Finite liquidity induces a higher sensitivity of the price 
increment (Walras' Law) to the excess demand. A study of a model 
for the order book would permit to relate the liquidity to the effective 
number of agents N. If one would know this dependence one could in- 
sert it in the coefficient of the Walras' Law, leading to an additional 
feedback of the dynamics with respect to N. 

As this example shows many other specific elements can now be ad- 
dressed in a systematic way. 

One point which is left out from our four basic ingredients is the ques- 
tion of the performance of the agents and their persistence in the mar- 
kets. Our agents can indeed change their startegies on the basis of the 
past price time series and on the herding phenomenon. On the other 
hand we limit the strategies to the two broad classes of fundamenlists 
and chartists. This is a drastic simplification with respects to models 
which permit a much broader choice of strategies which are selected by 
genetic or neural algorythms [4l[5l[6]. Our model shows, a posteriori, 
that the simplification in the two broad classes for the two strategies of 
stabilization and destabilization is enough to achieve a detailed under- 
standing of the SF and their self-organization. On the other hand, the 
possibility of a broader class of strategies and their choice on the basis 
of previous performance is certainly an interesting point. The problem, 
in this respect, is that if the price dynamics is highly simplified without 
all the details of an order book, the analysis of the performance can 
become highly misleading because unrealistic strategies could appear 
as very profitable. In our opinion, therefore, the question of the strat- 
egy selection on the basis of performance is an interesting one but it 
would require a much higher degree of realism in order to be properly 
addressed. 
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A How to define returns 

In economics returns at a A-scale are defined in general as log returns, 
that is 

r,+A=ln^ (45) 
Pt 
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instead in present papers (except for section!?]) we have defined returns 
as 

n+A = Pt+A - Pt- (46) 



In this appendix we discuss the reason for this choice. The log hy- 
pothesis is consistent with the idea that it is the logarithm of an asset 
price to perform a RW rather than the price itself. Morever in such a 
way the price is defined positive. In Fig. [26] we report the daily NYSE 
returns defined as difference of prices and as log returns. A simple 
visual inspection confirms the previous statement because only the log 
returns seem to be reasonably stationary. 
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Figure 26: Typical daily NYSE index. A visual comparison suggests that 
the muhiplicative hypothesis for prices seems to be better assumption (sta- 
tionary) rather than the linear one. 



Anyway we can always write the price at any time as pt+A ~ Pt + ^PA 
and, if the time scale A is quite small (from minutes to few hours), we 
can generally assume that ApA ^ Pt and consequently: 

r.,A = lnfl + ^U^^ApA (47) 
V Pt J Pt 

because pt is a slow variable at this time scale (see chapter 5 of |T5]). 
Equation (|47p suggests that the way to define returns in our model is 
irrelevant unless ApA is very large. For example in Fig. [27] we show 
the autocovariance function of the returns corresponding to the linear 
dynamics but considering different definiton of the returns themselves. 
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Figure 27: The autocovariance function of returns, of absolute returns and of 
squared returns do not depend on the defintion of the returns in our model 
because the two functions are nearly equal in both cases. In the insert 
we highlight that even in semi-logaritmic scale the functions are nearly the 
same. Here we consider the model with the linear dynamics. 



The results are extremely similar for the various definitions. 



B Samples and tails 

Let us suppose that the increments are independent and the increment 
pdf f(r) is symmetric and with zero mean (the assumption of inde- 
pendence is not realistic however it can give some insight). Let us also 
define the thershold \r*\ as the value above which we define an event 
as rare. The probability of a rare event is consequently: 

/ + 00 
/(r)dr :=/3(r*) (48) 

The probability Pr of observing at least one rare event among N events 
is: 

PriN, /?) = 1 - (1 - /3(r*))^ w 1 - e-^^('-*) (49) 

because generally /?(r*) <^ 1. This probabilty is approximately 1 for 
f3 > and it is nearly where /? ^ N~^. So (3{r*). fixes the value 
of N in order to observe rare events. In practice in order to observe rare 
events N should be larger than /3(r*)~^. So, if one considers relatively 
long times correlation the effective size of the sample (with respect to 
the time interval in question) becomes relatively small in terms of the 
number of effective steps of a RW. 
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C Pdf tails 



In Fig. [8] we have pointed out that the variance of returns does not 
enterely capture the scahng behavior of the whole return pdf. The 
motivations may be traced in the origin of the non gaussian tails due 
to the interaction between a pdf which tends to have constant width 
and one with increasing width. Fig. [5] suggest that tails are some- 
what streched by the gaussian with increasing variance (red curves). 
Conversely the central part, which gives the main contribution to the 
variance, tends to behave like the fundamentalist gaussian. Therefore, 
dividing returns with their variance makes that pdfs collapse into a 
unique curve in the central part. On the other hand, in order to have 
a collapse of the tails, the variance of the chartist gaussian would ap- 
pear as a more natural choice. On the other hand Fig. [28] gives only 
a partial confirmation of this picture. In fact the tails collapse in the 
same curve only for A < 1000 — 2000, wc recall that these non collaps- 
ing regime (i.e. large value of A) corrisponds to the region in which the 
approximation of section 5 is no more in good quantitative agreement 
with the simulated results. Therefore it seems that for small value of 
A, the tails grow as quick as the variance of the chartist case but for 
large A tails are slower than the chartist gaussian but faster than the 
fundamentalist gaussian. 




Figure 28: Here we show the pdf normalized with the variance of the case 
X = 1. The data collapse is not quite satisfactory showing that the scaling 
behavior of the non gaussian tails is related to the complex fluctuations in the 
population dynamics and it cannot be interpreted with simple arguments. 
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